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1 Introduction 

Recently progress on the understanding of regularization by dimensional reduc- 
tion has been achieved in three directions. First, a mathematically consistent def- 
inition avoiding the problem found in Refs. [1,2] was formulated, and a succinct 
method to check the symmetry properties of dimensional reduction was devel- 
oped [3], leading to the verification of supersymmetry in important cases at the 
two- loop level [4]. Second, explicit calculations demonstrated how dimensional 
reduction can be applied to multiloop calculations and how renormalization has 
to be carried out in a non-supersymmetric context [5] . This provides the basis of 
transition rules between various definitions of parameters such as Ug or m^ and 
is useful to derive the GUT-scale values of these parameters from the experimen- 
tal values [6] . Third, an obstacle in the application of dimensional reduction to 
hadronic processes was removed [7] by the resolution of the factorization problem 
of dimensional reduction found in Refs. [8-10]. 

The purpose of the present article is to further elaborate on the application of 
dimensional reduction to hadronic processes. In Ref. [7] we restricted ourselves 
to the case considered in Ref. [8], the real corrections to the process gg -^ ti, and 
showed that, despite first appearances, in the collinear limit these real corrections 
factorize into products of splitting functions and leading-order cross sections. 

Here we will consider real and virtual NLO QCD corrections to arbitrary 
hadronic 2 ^ (n — 2) processes with massless or massive partons. We will dis- 
cuss the infrared singularity structure and the associated regularization-scheme 
dependence of all these corrections, provide transition rules between the schemes 
and show that all singularities factorize. In this way we show that the framework 
of dimensional reduction is completely consistent with factorization, and we show 
how this scheme can be used to compute hadronic processes in practice. 

One of the main points of this article is the distinction of two different versions 
of dimensional reduction that have been used in the literature. One of the reasons 
why the factorization problem of Refs. [8-10] has remained unsolved for so long is 
that these two versions have mainly been applied by two different communities. 
The version used in Refs. [8-10] is the same as the one defined in Refs. [3, 11, 12] 
and is the one mainly used in the context of supersymmetry. The version used 
in Refs. [13-15], which was denoted by DR and is actually equivalent to the 
four-dimensional helicity (fdh) scheme [16] at one-loop, is mainly used in the 
context of QCD. For the latter version, the infrared singularity structure and 
transition rules have already been derived [13-15]. We denote these two versions 
by DRED and fdh. They differ in their treatment of external particles, in a way 
analogous to the difference of the "conventional" and " 't Hooft Veltman" versions 
of dimensional regularization, CDR and HV. 

In the main part of the present article we will provide results and transition 
rules for all these four regularization schemes, keeping in mind that the results 
for the FDH, CDR and HV schemes can already be found in Refs. [13-15], while 



the results for DRED are new. The results for the infrared singularity structure 
in CDR, HV, FDH and in DRED are developed in Section 2 and 3, respectively. 
The practical application of DRED and transition rules are discussed in Section 4. 
Appendix A provides explicit results for all relevant splitting functions, and in 
Appendix B we provide three explicit examples of NLO computations in DRED. 

1.1 Elements and scheme dependences of hadronic cross 
sections 

We consider n-parton processes with up to two hadrons in the initial state at next- 
to-leading (NLO) order in QCD. The partons can be either massless quarks q or 
gluons g or massive partons such as heavy quarks Q, gluinos, or squarks. In our 
equations we will restrict ourselves to the most interesting case of two initial-state 
partons as the simpler cases can be obtained by straightforward modifications. 
The cross sections of such processes can be written as 

da{Hi{Ki)H2{K2) ^ a3...an) = ^ dXifaj/Hi{xi) dX2fa2/H2{x2) 

... „^ Jo Jo 
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da[ai{xiKi)a2{x2K2); as . . .a„), (1) 



where Hi 2 are the initial-state hadrons, K12 their momenta and aj {i = 3 . . .n) 
the final-state partons. The sums run over all possible flavours of the initial-state 
partons 01^2 of the hard partonic cross section da, and fa./nX^) denote the appro- 
priate parton distribution functions (PDF). In the computation of hadronic cross 
sections, three scheme choices have to be made: the choice of the renormalization 
scheme, the factorization scheme, and the regularization scheme. 

The functional dependence of the hard cross section da on input parameters 
like as, particle masses, etc, depends on the renormalization and the factorization 
schemes. 

The choice of the renormalization scheme is equivalent to a precise definition 
of the input parameters entering the computation, in particular of as, particle 
masses and other coupling constants. Once it has been fixed and renormaliza- 
tion has been carried out accordingly, all off-shell Green functions are finite and 
unambiguously defined. Changing the renormalization scheme changes both the 
functional form of da and the numerical values of the input parameters, such 
that da is renormalization-scheme independent up to terms which are formally 
of higher order than NLO. Common renormalization schemes are the MS-scheme 
for as, the DR-scheme for super symmetric parameters, or the on-shell scheme for 
masses.^ In the following we are not concerned with the renormalization-scheme 
dependence and assume that some renormalization scheme has been fixed. 



"'^We stress that, although the MS- and DR-schemes have originally been defined with refer- 
ence to specific regularization schemes, they can be realized in the context of any regularization 
scheme. 



The choice of the factorization scheme is equivalent to a precise definition of 
the parton distribution functions. Both the functional form of the hard cross 
section and the numerical values of the parton distribution functions depend 
on this choice, but this dependence cancels in the full hadronic cross section 
up to terms which are formally of higher order than NLO. The major part of 
the following considerations is independent of the factorization scheme, and we 
simply assume that some factorization scheme has been fixed. At the end we 
will specialize to the important case of the MS-factorization scheme. In that 
case, the parton distribution functions fat/Hi can be taken for instance from the 
well-known MRST or CTEQ sets [17, 18]." 

Our main focus is the influence of the regularization scheme (rs). After re- 
moving the regularization, all quantities appearing in Eq. (1) are RS independent, 
but the hard partonic cross section is a sum of RS dependent parts. It is commonly 
written as 

da = da^r + dcr'^f + rf^^r + da^\ (2) 

where the RS dependence is explicitly indicated. The lowest-order, or Born cross 
section da^°"^ is finite, and in the limit where the regularization is removed its 
RS dependence vanishes. The three NLO contributions are the real and virtual 
corrections da^'^'^^, da^"^ and the collinear counterterm da'^°^^, which subtracts 
initial-state collinear singularities. All NLO contributions involve collinear and/or 
soft singularities and depend on the RS in their finite and their divergent parts. 
There are no ultraviolet singularities and associated RS dependences in da^^^^, 
because these are eliminated by renormalization and by fixing the renormalization 
scheme. 

1.2 Variants of dimensional regularization and dimen- 
sional reduction 

In all dimensional schemes space-time is continued from 4 to D dimensions, where 
D = 4 — 2e is an arbitrary complex number. In this way momentum integrals 
become well-defined and ultraviolet and infrared singularities appear as 1/e*^- 
poles as e — !► 0. Gluon fields are treated differently in dimensional regularization 
and dimensional reduction. In the former, gluons are treated as D-dimensional 
as well; in the latter, gluons are treated as 4-dimensional. 

Both choices have certain advantages. The purely D-dimensional treatment 
of all objects leads to simpler expressions, but it breaks supersymmetry owing to 
the different number of degrees of freedom of the gluon and the gluino. The 4- 
dimensional treatment of the gluon is better compatible with supersymmetry and 
it is more amenable to helicity methods, which are commonly used to simplify 
QCD higher-order computations. 

In order to formulate the two schemes one needs to distinguish three spaces: 



• the original 4- dimensional space (4S). 

• the formally D-dimensional space for momenta and momentum inte- 
grals. This space is actually an infinite-dimensional vector space with cer- 
tain D-dimensional properties [19,20], and is sometimes called "quasi- -D- 
dimensional space" (QDS). The space 4S is therefore a subspace of QDS. 

• the formally 4-dimensional space for e.g. gluons in dimensional reduction. 
This space has to be a superspace of QDS in order for the dimensionally 
reduced theory to be gauge invariant. Hence it cannot be identified with 
the original 4S — it can only be constructed as a "quasi-4-dimensional 
space" (Q4S) [3,21] with certain 4-dimensional properties. In practice the 
distinction between Q4S and 4S often does not matter, but it is important in 
the definition of the different versions of dimensional reduction, see below, 
and to avoid the inconsistency uncovered in Ref. [2]. 

These three spaces are characterized by their metric tensors, which we denote by 
gfiu ^^^^ Q4S), g'^'^ (for QDS), and g^'^ (for 4S). The dimensionalities of the spaces 
are expressed by the following equations: 

g^'g^.u = 4, g^'g^, = D = 4-2e, r^g.u = 4. (3) 

The following projection relations express that 4S is a subspace of QDS and QDS 
is a subspace of Q4S: 

g'^'g." = g"', g'''g/ = g'', g'^g/ = r'- (4) 

It is useful to introduce the orthogonal complement to QDS. This is a 4 — Z^ = 2e- 
dimensional space with metric tensor g^'^ , which satisfies 



g^'' = r''+r. 






(5) 


r^'g.u = A-D = 2e, 






(6) 


g^'g/ = r'. 


r^g/ = 0, 


r^g/ = 0. 


(7) 



Within this framework it is now possible to precisely state the calculational 
rules of dimensional regularization and dimensional reduction. Since momenta are 
always treated in D dimensions, it only needs to be specified how gluons (or other 
vector fields) are treated. More precisely, it needs to be specified which metric 
tensors are used in gluon propagator numerators and in gluon polarization sums. 
Particularly important for the understanding of factorization is the treatment 
of gluon polarization sums in squared matrix elements. Without regularization, 
such polarization sums can be written as 

pols 



internal gluon 
external gluon 



CDR HV FDH DRED 



j^fiu ^t^u ^t^u t^u 



Table 1: Treatment of internal and external gluons in the four different RS, i.e. 
prescription for which metric tensor is to be used in propagator numerators and 
polarization sums. For the definition of "internal" and "external" see text. 



where k is the gluon momentum and n is a gauge vector such that nk ^ 0. With 
regularization, the metric tensor in this polarization is replaced by either g^'^, 
gf", gf", or g^". 

It is not strictly necessary to regularize all gluons. Only gluons that appear 
inside a divergent loop or phase space integral ( "internal" ) need to be regularized; 
for all other gluons ("external") regularization is optional. The precise definitions 
of "internal/external" in this context are as follows: "Internal gluons" are defined 
as either virtual gluons that are part of a one-particle irreducible loop diagram or, 
for real correction diagrams, gluons in the initial or final state that are collinear 
or soft. "External gluons" are defined as all other gluons. 

Now, since external gluons do not have to be treated in the same way as inter- 
nal ones, it is in fact possible to distinguish two variants of each regularization. 
The two variants of dimensional regularization are: 

• CDR ("conventional dimensional regularization"): Here internal and exter- 
nal gluons (and other vector fields) are all treated as D-dimensional. 

• HV ("'t Hooft Veltman scheme"): Internal gluons are treated as D- 
dimensional but external ones are treated as strictly 4-dimensional. 

Note that the above definition of internal gluons in phase space integrals is nec- 
essary for unitarity but leads to complications in the treatment of phase space 
integrals in schemes where internal and external gluons are treated differently. 
The two analogous variants of dimensional reduction are: 

• DRED ("original/old dimensional reduction"): Internal and external gluons 
are all treated as quasi-4-dimensional. 

• FDH ("four-dimensional helicity scheme"): Internal gluons are treated as 
quasi-4-dimensional but external ones are treated as strictly 4-dimensional. 

Table 1 illustrates these four schemes. 

Note that the version of dimensional reduction denoted by DR e.g. in Refs. [13, 
15] is equivalent to FDH at the one-loop level (see e.g. Refs. [13,22]).^ The infrared 



^In Rcf. [23] a two-loop definition of the fdh scheme has been given. In what follows we 
will only use the one-loop definition and the one-loop equivalence of fdh and dr. 



properties of the three schemes CDR, HV, DR (or equivalently fdh) have been 
studied and compared in Ref. [15] and found to be consistent with factorization. 
An apparent inconsistency between dimensional reduction and factorization 
has been identified in Refs. [8,10], but in these references the version dred has 
been used. In Ref. [7] it was found that factorization holds as expected in dred 
if external quasi-4-dimensional gluons are decomposed into D dimensional gauge 
fields and (4 — D) dimensional "e-scalars", which are treated as separate partons. 
Technically, this decomposition amounts to replacing 

•Mdred(- ■■g---)= -^dred(- • • ^ • • •) + -^dred(- • • ^ • • •) 

= E -^dr,ed(---^---), (9) 

ffG{g,3} 

for squared matrix elements, where the different gluon types g^ g, g are denoted 
by the same symbols as the associated metric tensors. The algebraic expressions 
for the partonic processes involving g, g, or g are defined by the values of the 
corresponding gluon polarization sums. These read 



pols 
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^ (nk) {nkf ' ^ ' 

~g, ^e^e-^-^^^ (10c) 

pols 

Eq. (9) follows trivially from Eq. (5). 

The decomposition of gluons into their D-dimensional and e-scalar part in 
DRED is also relevant for the renormalization of UV divergences. In order to make 
Green functions with external e-scalars finite, the renormalization constants for 
e-scalar couplings in general have to be different from the corresponding gluon 
couplings. For example, in pure QCD, the couplings a^ and «£ for the quark- 
antiquark-gluon and the quark-antiquark-e-scalar vertices receive different coun- 
terterms 5as 7^ 5ae even if as = fte at tree level [5,24]. 

1.3 Splittings in the four schemes 

An essential part of the RS dependence of NLO contributions is related to the 
RS dependence of the splittings i -^ jk of one parton i into two collinear partons 
j, k. The RS dependence of real corrections is related to the splitting functions 
^i^jk'i the RS dependence of virtual corrections is related to constants 7rs(0 [13], 
which in turn can be derived from the P^jj^ via unitarity [15]. In this section we 
explain the RS dependence of the splitting functions and correspondingly of the 







CDR 



HV 



FDH 



DRED 



Figure 1: Gluon splitting into two collinear gluons in the four schemes, indicating 
the appropriate treatment of each gluon. 



7Rs(i). The full results can be found in Appendix A. Figure 1 shows the most 
interesting case of a gluon splitting into two collinear gluons. According to the 
definition given above, the two collinear gluons j and k are treated as "internal" , 
and the virtual gluon i as "external" . The appropriate treatment of the gluons 
in the four RS can be read off from Table 1 and is displayed in the figure. Two 
simple observations allow an easy comparison of the four cases. 

First, the projection of a D-dimensional onto a strictly 4-dimensional parent 
gluon does not change the structure of the result of the splitting functions. And 
second, the result in DRED should be decomposed according to Eq. (9) into four 
splittings g — > gg, g -^ gg, g -^ gg, g -^ gg.^ Then the result in CDR is identical 
to the DRED result for g -^ gg, and all scheme differences can be explained in the 
following way: 

• The splitting g ^ gg is identical in the CDR and HV schemes. Formally, 
this is expressed in the equality 



p<CDR/ X 



P 



<HV 



■*99^ 



p<DREOf \ 



;ii) 



for the splitting functions defined for z < 1. 



In the FDH scheme the outgoing gluons are treated as quasi-4-dimensional. 
The resulting additional term can be interpreted as being due to the split- 
ting g ^ gg as already discussed in Ref. [15]. Hence, 



P^Z7M = Pg^ZTi^) + p^, 



)<DRED|' 



(12) 



In the DRED scheme the parent gluon is also treated as quasi-4-dimensional, 
and therefore the two additional splittings g ^ gg and g ^ gg are possible. 
In the spirit of our discussion around Eq. (9) we do not combine the DRED 
splitting functions into a single one. 



^Splittings involving an odd number of g vanish. 



The splitting functions involving quarks are related in a similar way. Via uni- 
tarity, the RS dependence of the constants 7(2) follows from the splitting func- 
tions [15] and can thus be explained in an analogous way: 

• The 7(2) in CDR and HV are the same, 

7cdr(0 = 7hv(0 for ^ e {9, q}- (13) 

• The additional terms in the fdh scheme are due to the splittings g ^ gg 
and q -^ qg: 

7pdh(^) = 7hv(^) - fdzzP^^l^-{z), (14) 

7POH(g) = 7Hv(g) - fdzz [P<^^r{z) + P<4r(^)] • (15) 

Jo 

• In DRED one has to distinguish 7 constants for g, g, and q. The ones 
corresponding to g and q are the same as the ones in the FDH scheme; the 
one for g is related to the additional splittings g -^ gg and g — * gg: 

7dREd(^) = 7FDH(fi'), (16) 

7DRED(g) = 7FDH(g), (17) 



7dred(^) = - / dzz- — [P-g 

Jo U ^ ^)+ 



<DRED/ \ 



+Pg^^ri^) + 2iv^ p^4r (^)j (18) 

These relations form the basis for understanding the RS dependence of NLO con- 
tributions and in particular the difference between DRED and the other schemes. 
In the subsequent sections we will see that additional RS dependences arise from 
the crossing of the splitting functions to initial-state parton splitting and from 
the RS dependence of the LO matrix element. 

2 CDR, HV, FDH 

Our starting point is the decomposition, Eq. (2), of the hard partonic cross sec- 
tion, and we are mainly interested in the RS dependence of the separate terms 
contributing to da. In this section we will restrict ourselves to the well-known 
cases of CDR, HV and FDH. The RS dependence of quantities will be indicated 
by a subscript rs*, the star reminding us that we consider CDR, HV and FDH, but 
not (yet) DRED. 



2.1 Born term 

We consider the partonic process 



ai(pi) 02(^2) -^ asiPs) ■ ■ ■ an{pn) 



(19) 



where a^ and pi denote the flavour and the momentum of parton i respectively. 
The /-loop correction to the RS dependent squared matrix element for the pro- 
cess given in Eq. (19) is denoted by A<Rs*(ai(pi), 02(^2); ^3(^3) • • -CiniPn)) or by 
A1rs*(«i • • ■O'n) for short. For the cross section we need the averaged squared 
matrix elements 



til) 



(Aiy*(ai,a2;...a„)) 



1 



1 



,(0 



2 S12 u;Rs*(fli) <^Rs*(fl2) 



A<y*(ai,a2;...a„), (20) 



where cJrs* (aj) denotes the RS dependent number of degrees of freedom of a parton 
with flavour a-i and S12 = 2 (pi ■ P2) in the case of massless incoming partons. 

The Born cross section is obtained by integrating the squared and 
averaged tree-level matrix element over the (n — 2) parton phase space 
(i$n-2(pi,P2;P3 • • -Pn) multiplied by a measurement function for an infrared-safe 
quantity and a symmetry factor. The latter two are always implicitly understood 
in our notation and we simply write 



(iaRsl(ai . . .a„ 



/ d^n-2{pi ■■■Pn) (-MRsl(ai • • • a„)). 



(21) 



The RS dependence in Eq. (21) is due to 0{e) terms in A^rs*. Since we consider 
an infrared-finite quantity, the phase-space integration does not introduce any 
poles. Therefore, we can take the limit e — > and 



da 



Born/ 



(ai . . . a„) = da^^J,^{ai . . . an) 



_D->4 



(22) 



is RS independent, as indicated by the absence of the subscript rs*. 



2.2 Virtual corrections 

For the virtual corrections we need A^rs*, the interference terms of the one-loop 
amplitude and the tree-level amplitude. The structure of the singular terms of 
AIrs* is well known [14,25]. For the fully renormalized matrix element, it is given 
by 



M^nL{ai ■ ■ ■ ttr, 



2^ 



cr 



,(0) 



MUMl ■ ■ ■ ^r, 



( --Xl^i^s*(ai)] 



J2 V(«, j) M%,{ai . . . a„) + M'i^^ (ai . . . 



«j 



(23) 



where the sums i and j are over all initial or final state partons and we introduced 

The soft and coUinear poles are contained in the terms proportional to A^rs* and 
A^Rs*- The latter are the colour-linked Born squared matrix elements introduced 
in Ref. [26] and correspond to the square of the colour-correlated tree amplitudes 
with a Tj ■ Tj insertion, used in Ref. [27]. If particles i and j are massless, we 
have 

V(,,).__LRe(_^)-' ,25) 

and Eq. (23) reduces to the well-known expression for the singularities of one- 
loop QCD amplitudes [25]. If one or both of the particles i, j are massive, these 
expressions have to be generalized [14], but the structure of the singularities 
remains as in Eq. (23). 

The RS dependence of Eq. (23) is contained in the constants 7Rs*(aj)/e, as 
well as in M-H* and TWrs*. The remaining term, A^ns(oi . . .On), is in general 
very complicated, but is finite and, after taking the limit -D -^ 4, RS independent. 

The essential, non-trivial part of the RS dependence is due to the 7Rs*(ai)/e 
terms. These are closely related to collinear singularities due to self-energy inser- 
tions on external legs and depend on the fiavour a, of leg i. The RS dependence 
of 7rs* has been given in Ref. [13]. Via unitarity it is related to parton splittings 
tti — i^anything [15], as summarized in Section 1.3. 

In the present paper we determine 7rs* by insisting that the sum rules hold in 
all RS to all orders in e. This is a slightly different approach compared to Ref. [13] 
and simply amounts to a shift of finite terms between Aiiii and the 7rs* terms 
in Eq. (23). Neglecting O(e^) terms, we find 

, , , , /3o , TpNp , .. , .. SCp , Cf 

1cdbX9) = 1hy{9) = y + ^ g ; 7cDR(g) = 7Hv(g) = ^- + e — 

/?o ^ 2TfNf - N, 3Cf 

7pdh(^) = y + ^ Q 5 7FDH(g) = ^- (26) 

with Po = (llA^c — 4:Tf Nf)/S and Tf = 1/2. For heavy quarks the result is 
RS independent [14] and we have jiQ) = Cf- 

The final virtual corrections da^l are obtained as 

da^l^liai . . . a„) = / d^n-2iPi ■■■Pn) {M^nliai . . . a„)). (27) 

Since this phase space integration does not give rise to any pole in e, taking the 
hmit D ^ 4 in Aij^i is justified. 
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2.3 Real corrections 



For the real corrections da^^g^: to the partonic process ai a2 -^ as . . . a„ we have 
to consider contributions from all 2 — > (n — 1) processes that are obtained by 
a split of any of the outgoing partons. We have to evaluate the corresponding 
squared matrix elements A^^g (ai, 02; 03 . . . dn+i) and integrate them over the (n — 
1) parton phase space d^n-i{Pi,P2'iP3 ■ ■ -Pn+i) 

^C^rT*' = 5Z / d^n-l{Pl,P2;P3 ■ ..Pn+l){M^RS*{ai, ^2] ^3 • • • On+l))- (28) 

ai 

In Eq. (28) we denote by Oj, i G {3 . . . n + 1} the flavour of the outgoing partons, 
and as indicated by X^a ' ^^ have to sum over all relevant processes. 

As is well known, the matrix elements can develop singularities in regions of 
the phase space where a parton becomes soft or two partons become collinear. 
The integration over d^n~i{pi ■ ■ -Pn+i) in this region then results in 1/e^ and 1/e 
poles. Thus the RS dependence of the matrix elements which manifests itself in 
the 0{e) terms of A^rs results in differences in the 0{l/e) and in the finite terms 
of the real corrections.^ 

In order to deal with the phase-space integration at NLO one often uses either 
phase-space slicing [28] or subtraction [26,29], and several general procedures have 
been developed [27,30,31]. For our purposes it is sufficient to know that they all 
rely on the same main points. Using the results given below it will be obvious 
how any of these procedures can be applied in the context of dred. 

The first point is that in any of the singular regions the matrix elements take 
a simple form and can be written as a factor containing the kinematic singu- 
larity times a reduced (colour-linked) tree-level matrix element, associated with 
a 2 ^ (n — 2) process. Secondly, the phase space is factorized according to 
d^n-i{Pi ■ ■ -Pn+i) = d^n-2{Pi ■ ■ ■ p'n) '^'^rad- The factor with the kinematic singu- 
larity is integrated analytically over d^j-ad, producing the poles in analytic form. 
These poles will be multiplied by the reduced (colour-linked) matrix element 
and are to be integrated over a slightly modified {n — 2) parton phase space 
d^n-2{Pi ■ ■ -Pn)- It is therefore not surprising that the real corrections have a 
similar structure as the virtual corrections, Eqs. (23) and (27). We will now look 
at all three potentially singular regions in turn. 

Soft Region: In the limit where gluon g^ (or another massless gauge boson) 
becomes soft we have 

M'aliai, as; . . . gk{Pk) ■ ■ ■ a„+i) ^''= g] ^ — '■^— A^r^s*(«i •••««)> (29) 

• • ^ikSjk 



''Note that in the HV and fdh schemes real soft and/or cohinear gluons have to be treated 
as "internal", i.e. in the same way as gluons in a closed loop but differently from observed, 
"external" gluons. This is the only source of the RS dependence of the real corrections in these 
schemes. 

11 



where we introduced g1 = Ait as- It is understood that the set of partons 
{03 . . . a„} is equal to the set {03 . . . g^ . . . a„+i} with g^ removed. Similar com- 
ments apply to the analogous equations below. The phase space integration of 
Eq. (29) leads to J d^ra.dSij/{sikSjk) and results in poles that cancel the corre- 
sponding poles in V{i,j). Due to the measurement function implicitly included 
in d^n-2, the remaining integration does not result in any singularities. The 
scheme dependence enters only through A^rs* and after summation over all real 
processes trivially cancels between the real and virtual corrections. 

Final-State Collinear Region: In the limit where two outgoing partons dk and 
0,1 become collinear we have 

M'aLiai, as; . . . ai{pi) . . . ak{pk) . . . a„+i) ^=^' (30) 

9 rfl 

■^ ys p<RS* (^\ »^(0) 



Ski 



■P(k^**^ki(^) -^Rs*(ai, 02; . . . a(M)(pfc + Pi)--- an). 



As detailed in Appendix A, -P^^^^^^^ is the RS dependent splitting function defined 
for z < 1 with pk — > z{pk +Pi) and p; ^ (1 — z){pk +Pi)- In Eq. (30) the flavours 
Ofc and di are flxed. This uniquely determines the flavour of the parent parton 
a(fci). It is understood that if the split is flavour forbidden we set -P/^n'^!^^; = 0. 
To avoid a proliferation of subscripts, we denote the flavour of the partons in the 
splitting functions simply by {kl) etc. rather than a(^ki)- The parent parton is 
slightly off shell as indicated by the notation (kl)*. 

Contrary to the soft limit, in the collinear limit there are two sources of 
RS dependence. Apart from the trivial dependence through A^Rg*(ai . . . a„), the 
prefactor -P/^n.^^^^ is also RS dependent. Its RS dependence can be found in 
Section 1.3 and in Appendix A. Since the z dependence in Eq. (30) is entirely in 
the prefactor, the integration f rf'I'rad P(ki)*l^ki/^ki can be performed separately. 
The terms related to the collinear singularities due to the splitting of parton 
o-i = a(ki) schematically can be written as 

d(y'^^}\ai . . . ai . . . an) = --^— rfcri°\(ai . . . a^ . . . a„) (31) 

ZTT e 

where the sum over all possible splittings, Yla ' ^^ ^^^ to the sum over the relevant 
real processes, Eq. (28). Note that in the sum Ofc G {g,q,q} a sum over the Np 
massless quark flavours is implicitly understood. After this sum the integrand is 
symmetric with respect to z ^^ 1 — z. Hence the integration can be restricted to 
the region z > 1/2 and the potential singularity at ^ = 1 is regularized with the 
usual +prescription. The z integration in Eq. (31) results in a RS dependence of 
the singular prefactor that multiplies the Born term. In fact, the factor in the 
second line of Eq. (31) is equal to (— 7Rs*(aj)). Therefore, after summing up the 
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7rs*(^)L„ ^ _ 7rs*(c/)|^^ 



+ 





Figure 2: Illustration of Eq. (31) for the case where Oj is an outgoing splitting 
gluon. The sum over all relevant real processes, ^^ , gives rise to two contribu- 
tions. The one on the left (right) results in the N,. [Np) part of 7Rs*(fl')- 



contributions of all final state partons, i G {3 . . . n} these terms precisely cancel 
the singularity and the RS dependence of those virtual terms displayed in the first 
line of Eq. (23) that are associated with outgoing partons, i > 3 [15]. 

Initial-State Collinear Region: Finally we turn to the case of an outgoing par- 
ton afc becoming collinear to the incoming parton ai (or 02). There are some im- 
portant differences with respect to Eq. (30). To start with, the collinear limit has 
to be written with the spin/colour summed/averaged matrix elements (A^rs*). 
With pfc -^ (1 — z)pi the collinear limit is given by 

{Mns*{ai{pi), 02; . . . afc(pfe) . . . a„+i)) ^=^' (32) 

Contrary to Eq. (30), the z dependence in Eq. (32) is not restricted to the pref- 
actor and, therefore, the ^-integration results in a more complicated structure. 
In Eq. (32) the flavours ai and dk are fixed and uniquely determine the flavour of 
parton a^ik), which is slightly off shell. The splitting functions of Eqs. (30) and 
(32), with initial/final-state off-shell parton, are related by a crossing relation 

p,<4K,..(^)-(-i)*-^^-^S-.(i). (33) 

where #/ denotes the number of crossed fermions. We remark that the well- 

CDR p< CDR 

*{lk)*k ~ ^l*^(lk)k 



known crossing symmetry P^n^,i^ = Pif^^ik)^ does not hold in all RS, and in 



Pitl^vki^) = P^f^'mki^) + ^r:mkiz) (34) 



general we have 

^ *{ik)*ky^) = ^i*^{ik)kK'^) "^ ^i^{ik)k\ 

with A™^ ^ and A"^ =^ 0. The explicit form of A''^" and A"^ can easily be 
found using Eq. (33) and the results in Appendix A. However, as will be discussed 
in Section 2.4, the distinction made in Eq. (34) is ultimately not required. 

Using Eqs. (34) and (86) to express the collinear limit, Eq. (32), in terms of 
the full splitting functions -Pfil*^^-)^ and summing over all relevant real processes, 
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Figure 3: Illustration of Eq. (35) for an incoming splitting gluon. The sum over 
all relevant real processes, Y2a ' gives rise to three terms. 



we can write the initial state collinear term for parton 1 schematically as 



^^RsV (ai-.-an) 



a^cr 
27r e 



7rs*('^i) dans*{ai{pi), 02; • ■ ■ On) 



(35) 



The first term on the r.h.s. of Eq. (35) is due to the 6{1 — z) term present in 
Eq. (86). Together with the corresponding term for the second incoming parton, 
da^^s* ) this results in a term that precisely cancels the singularity and the RS 
dependence of those virtual terms displayed in the first line of Eq. (23) that are 
associated with incoming partons i <2. The remaining terms given in the second 
line on the r.h.s. of Eq. (35) are associated with collinear counterterms. 

2.4 Collinear counterterm 

In the sum of the virtual and real corrections, da^ll + da^^^l^ all singularities and 
RS dependences cancel, apart from the terms given in the last line of Eq. (35). 
These are cancelled by the collinear counterterm da^^l- While the divergent parts 
of the collinear counterterms are completely determined, there is some freedom in 
how to specify the finite parts of da'^^l- Any specific choice of the finite parts of 
da'^^l is equivalent to the definition of a particular factorization scheme. Leaving 
the factorization scheme open, we can write 



'^^RS*,Fs('^l' '^2; • • • «n) — -^ 2^ / dz 



(36) 



X 



(P-*,(^) + A«- ,(^) + eXll,,{z)) da'^l{a,{zp,), a,{p,)- . . . a„) 



(0) 



(0) 



+ {P2%k{z) + A«!:,fc(z) + tXll^^{z)) da'^l{ai{pi), ai{zp2)] . . . a„) 



where the sum is over all possible splittings of the incoming partons, and the 
index i is defined in analogy to Eqs. (30) and (31). The Xfj^-/^ are the finite 
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(i.e. e-independent), RS independent terms which define the factorization scheme. 
The formulas in Sections 2.2 and 2.3 show that in this way the hard (subtracted) 
partonic cross section 



ddpsitti . . . a^ 



rfCr(ai...a„) (37) 



+ dcr™l{ai . . . a„) + da'j^g^{ai . . . a„) + c?o-^s*,fs('^i • • • «n) 



D^4 



is finite and RS independent as indicated by the absence of the subscript as. In 
this approach the most natural factorization scheme would correspond to setting 
all X^^ = 0. In principle, this particular factorization scheme is as good as 
any other, and it could be used in practice in Eq. (1) in conjunction with parton 
distribution functions fa,/H determined in the same scheme. In practice, however, 
parton distribution functions such as the standard MRST or CTEQ sets [17, 18] 
are mainly available in the MS factorization scheme, which is different. 

The MS factorization scheme is defined by using rs=CDR and replacing the 
square bracket in Eq. (36) by 



[PS-fc(^)]z5^4 c^41(ai(^pi), aa; . . . a„) + {1 ^ 2} 



(3^ 



The minimal subtraction procedure corresponds to setting e ^ in the splitting 
functions PiZfk- Thus, even if we use CDR (and even though A"^^/^ = 0) the MS 
scheme does not correspond to X^^ = but to 

eX^,,{z) = ~P^!,{z), (39) 

where P^-,,{z) = Pf_^f;,(z) - [Pi^fkiz)]^^^^^ denote the 0(e) terms of the splitting 
functions in CDR. 

In the evaluation of the real corrections, Eq. (35), as well as in the collinear 
counterterm, Eq. (36), we made the distinction between -P;f_!^ffc and PiIl^^I, the 
splitting functions appropriate for an outgoing and incoming split, respectively. 
In a general RS, these two splitting functions differ as indicated in Eq. (34), hence 
the presence of the A terms in Eqs. (35) and (36). While the expressions given in 
Eqs. (35) and (36) are those that naturally arise in the calculation, we note that 
the A terms cancel in the sum of da^j^^^ + da^^lpg. Thus we can drop A^^^^-^^ 

in Eq. (35) (and the corresponding term in da^g^' ) if we also drop A^^-/^ (and 
A^!^*.;.) in Eq. (36). Ultimately, the distinction between P^fJ^f^ and P^l^f^ is not 
needed. 



3 Dimensional Reduction 

In this section we show how the structure described in Section 2 can be generalized 
to include dred. As mentioned in the introduction, the key point is to split the 
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gluon into a D-dimensional gluon g and a (4— D) -dimensional e-scalar g by setting 
9 = 9 + 9 and using Eq. (9). Often it is sufficient to perform this split for only 
one gluon as in Eq. (9); in general, if a process contains 4^g gluons, gi^ . . -gi^g, 
we can decompose the matrix element in dred into 2*^ terms according to 

^DREoi- ■ ■ 9h ■ ■ ■ 9i#, ■ ■ ■) = Yl ■■■ 5^ '^DRED(---^n---^i#3---)- (40) 

We consider the two partons g and g to be two different partons and consequently 
regard the r.h.s of Eq. (40) as a sum over the squared matrix elements of 2*^ 
different processes. To bring our notation in line with the previous section, we 
will write Eq. (9) and Eq. (40) as 

•^dredIoi •••«*•••«") = 5]]A^DRED(oi---ai...a„), (41) 

-^DRED(«l---«n) = ^A<DRED(ai---an), (42) 

{»} 

respectively, where it is understood that li ai = g we sum over the two terms 
^j ^ {^) 9} whereas if Oj = g there is only one term in the sum Sj G {q}. 

For the spin summed/averaged matrix elements the relation equivalent to 
Eq. (42) reads 

/.. / x\ V^ ^DRED(fll) ^DRED(a.2) /■ . /v v v n\ /<o\ 

TTT ^DREdIOiJ '^DREDla2J 

The explicit expressions for ui-nsio-i) are given in the Appendix in Eq. (85). 

We stress that while the split g = g + g is conceptually simple, it seems to 
complicate practical computations. As we will see in the later sections, however, 
in an explicit computation of a physical process in dred it is only required at a 
very limited number of steps. In particular, it will turn out that no PDF for the 
unphysical e-scalar g will be required. 

In the present section we will use the split to understand the infrared structure 
of matrix elements in DRED. In fact, it is straightforward to see that Eqs. (23), 
(30), and (32) for the collinear singularities of virtual and real corrections hold in 
the same form in DRED for the individual processes with split partons, i.e. if we 
replace rs* -^ dred and Oj -^ hi in these equations. However, our main interest 
are the infrared properties and RS dependences of matrix elements for full gluons; 
therefore we will carry out the sums over {cii} wherever possible. 
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3.1 Born term 

The full tree-level matrix element in dred is equal to the one in fdh and HV and 
can be obtained from the CDR result simply by setting D ^ 4 

-MdredIch ■■■an) = Mnliai • • • «„) = M'p^n{ai ■■■an) (44) 

= M'-c^niai ■ ■ ■ an) 

L J 1)^4 

The Born cross section in DRED can be obtained from Eq. (44) and satisfies 

da^uREoiai ...an) = da^^viai ...an)= daiuniai ...an)^ (45) 

3.2 Virtual corrections 

The structure of the virtual corrections in DRED is analogous to Eq. (23) for 
each individual DRED process with split partons Oj. Hence, by summing over all 
processes as in Eq. (42) we obtain 

M'hLniai ■■■an) = y^ TT^cr 

{a} 



-MDRED(ai...a„) f — ^ 
+ J2 ^(^' •^■) -^dred(«1 •••««)+ M^'i (ai . . . 




The non-trivial structure of Eq. (46) is most essential for the 7 terms, since 
7dred(^) 7^ 7dred(^)- In fact, as discussed in Section 1.3 and Appendix A, the 7 
for g and q match the ones in the FDH scheme, while the one for g is new and 
different: 

7dred(^) = 7fdh(^) = Y + ^ Q ' (47a) 

7DRED(g) = 7FDH(g) = ^^, (47b) 

7dred(^) = 2N, - TfNf. (47c) 



On the other hand, V(i,j) is not affected if a ^ is replaced by a g, and the 
Born terms in the last line on the r.h.s. of Eq. (46) can be combined in a straight- 
forward way. Hence one can immediately obtain the result for either the fully 
or the partially combined process, A^DRED(ai • • • ^n) or AlDRED(fli, ^2; 03 . . . a„), if 
desired. In particular, the DRED result for the process involving only full gluons 
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satisfies 



AlDRED(ai...a„) = --cr 

ZTT 



y^ y^ -A/iDRED(ai . . . Sj . . . a„) ( — TdredIoj) j 

(4J 



+ X] V(i, j) A<d^red(«1 • • • «n) + A<NS (Ol . . . an 

where A^dredIoi • • • On) = -^FDH(ai • • • ««) and A^y (ai . . . a„) is the RS indepen- 
dent term appearing also in Eq. (23). 

We recall that A^dred denotes the fully renormalized one-loop matrix element 
and Eq. (46) does not contain any ultraviolet singularities. In DRED this implies 
that off-shell Green functions are finite also if external g are present. This requires 
that couplings involving g and couplings involving g in general renormalize dif- 
ferently [24]. As a result, the renormalization procedure in non-supersymmetric 
theories can be slightly more involved in DRED. We refer to Appendix B for 
examples. 

3.3 Real corrections 

The calculation of the real corrections in DRED follows the same pattern as in 
the other schemes discussed in Section 2.3. However, there are some important 
differences which we will consider for the three singular regions in turn. 

Soft Region: In close analogy to Eq. (29), we have to consider the limit of 
the real matrix element when a gluon becomes soft. In DRED this soft gluon 
can be either a, g oi & g. The soft limit is governed by eikonal factors of the 
form p'^p'^/{sikSjk), contracted with the corresponding polarization sum of the 
soft gluon gk- The polarization sums in Eq. (10) then show that a soft cjk leads 
to the same limit as a soft full gluon gk, while a soft e-scalar gk leads to zero. 
Hence, 



A^DRED(ai, as; . . . gk{Pk) ■ ■ ■ a„+i) *""= fi-^ ^ — '-^ -^dred(«i • • • «n), (49) 
M^uLuiai, as; . . . hiPk) ■ ■ ■ On+i) ^'= 0. (50) 



In analogy to the behaviour of V{i,j) in the virtual corrections, the soft limit 
Eq. (49) does not require the split gk = Qk + Qk- Thus, with respect to the soft 
limit, DRED is equivalent to HV and fdh. 

Final-State Collinear Region: Here the split g = g + g is essential. According 
to the main result of Ref. [7], the key equation for the collinear limit, Eq. (30), 
has to be modified in DRED if the parent parton a(^ki) is a gluon. In this case the 



Tdred 




Figure 4: Illustration of Eq. (52) for the case of an outgoing splitting gluon. The 
sum over all relevant real processes, ^^ , together with the sum due to the split 
Oj G {9^9} gives rise to four terms, resulting in the Nc and Np parts of 7dh,ed(5) 
and 7dred(5) respectively. Gluons g and (anti)quarks are drawn as usual. Dashed 
lines represent g and g is represented by a zigzag line. 



flavour of the parent parton is not uniquely determined by a^ and ai and we have 
to sum over the two possibilities a^ki) £ {g,g}- Thus, in dred Eq. (30) becomes 

MoLuiai, 02; . . . ai{pi) . . . ak{pk) ■ ■ ■ a„+i) ^=^' (51) 



2 0^ 

Yl ^(H™ (-2) -^DRED(ai, 02; . . . a^ki){Pk + Pi).-- an). 



«{fci) 



We remark that Eq. (51) leads us to consider splitting functions with a(^ki) = 9 
and a(fci) = ^ in DRED. They are given in Appendix A. Eq. (51) remains true if 
all partons aj are replaced by a^, however, only for a^ = 0,(^^1^ the split is strictly 
required. Thus, Eq. (31) is modified to 

da^j^HEoif^i ■ ■ ■ (^i ■ ■ ■ (^n) = --^— y^da'^lT,j,{ai...ai...an) (52) 






This is illustrated in Figure 4, where we consider the particularly interesting case 
of an outgoing splitting gluon. For each choice Oj G {g, g} we have to sum over 
all splittings a^ G {g,q,q}- Thus, the two terms of Figure 2 become the four 
terms of Figure 4. 

Initial- State Collinear Region: As for the final-state coUinear region, the key 
point is the necessary split g = g + g in the factorization of the collinear limit. 
In DRED, Eq. (32) has to be generalized to 

(A^DRED(ai(pi), 02; • • • dk{pk) ■ ■ ■ dn+i)) ''=^' (53) 

^ J2 ^il°S"^(^) (M^uLnidmi^Pi)^ «2; . . . a.)). 
S(ife) 
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Figure 5: Illustration of Eq. (54) for an incoming splitting gluon with partons 
represented as in Figure 4. The sum over all relevant real processes, Y2a ' together 
with the sum due to the split a(ijt) e {9i9} gives rise to four terms. 



Again, as far as the collinear limit is concerned, not all gluons in Eq. (53) have 
to be split. Only for the virtual parton a(ifc) the split is essential. 

Note that in dred the crossed splitting functions satisfy the crossing relation 
Eq. (34) without A terms. Hence, the initial state collinear term for parton 1 
can be written as 



7 real, 1/ ^ \ _ Ois Cr 

ZTT e 



y^7DRED(Qi) daoI.ED{aiipi) , 02; . . . a^) (54) 

0.1 

- ^^ d,z Pn™fc)fc da^oLD{a^ik){zpi), Oa; . . . a„) 

As illustrated in Figure 5, in the sum over 0(1^) in Eq. (54) it is essential that we 
treat g and g as separate partons, whereas a^ G {g, q, q}. For an incoming gluon, 
the three terms of Figure 3 are generalized in dred to the four terms of Figure 5. 

3.4 Collinear counterterm 

The collinear counterterm in DRED can now be constructed in the same way as 
in the other schemes. Generalizing Eq. (36) to DRED we can write 

c^^drUfs(«1) «2; • • • an) = ~^ J2J2 h^^ ^^^'> 

Ok Oi 

X (Pf^™(^) + ^Xll^^{z)) da^^Bnid^izpi), 02(^2); ...an) 
+ {P2^1k{z) + eX^^ii^iz)) rfcr^RED(ai(pi), di{zP2)] ...a„ 

The Xl^j^j^ involving partons Oj G {g, g} appearing here are defined in terms of 
the functions appearing in Eq. (36) as 



UJr 



•yi-'S — iJixi^D\g ) x^FS 



UJr 



ig) 



(56) 
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With these definitions, the resuhs of this section, Eqs. (48) and (54), show that 
the hard partonic cross section 



dapgitti . . . a^ 



rfC^,'^(ai...a„) (57) 



+ iio"^'^j,n(ai . . . a„) + da^^-i^^{ai . . . an) + c?o"dred fsI*^! • • • ^n) 



D^A 



is equal to the one in the other schemes given in Eq. (37). This shows in particular 
that it is possible to realize the MS factorization scheme in dred in the same way 
as in CDR, HV, or the fdh scheme. In order to make this result explicit we close 
the section by providing the full form of the appropriate collinear counterterm, 
valid in all RS, 

C^<°s"Ms(«b «2; • • • On) = ^ Y 5Z 5Z / ^^ (^^) 

{P^Z^ki^) - Pf^tkiz)) da^XUzpi).a,- . . .) + {1 - 2} 

where the sum ^^^ runs over ak G {(?, g, q} in all RS, whereas the sum ^- . runs 
over in di G {(/, g, g} in CDR,HV, FDH, and over di G {^, g, q, q} in dred. Also, in 
DRED we define^ 



jMSe /^\ — ^dREd(^) 



Pi^Mki^) = -^^^^^ PiZU^) (59) 

in analogy to Eq. (56). 



4 Applying DRED 

In the previous two sections we discussed how the singularities and RS dependence 
between the various parts of Eq. (2) cancel, and we found that the subtracted 
partonic cross sections given in Eqs. (37) and (57) are finite and RS independent. 
This is precisely what we wanted to achieve. However, there is still one conceptual 
issue to be addressed. 

The question is whether in the convolution of subtracted partonic cross sec- 
tions with PDF we need to distinguish between g and g in DRED. We will show 
that this is not the case. This will also entail that no PDF for finding an unphys- 
ical ^ in a hadron will be required. 

Once this issue is clarified, we will summarize our results and give transition 
rules between the various RS separately for all parts of the subtracted finite 
partonic cross sections. 



^The A terms appearing in the HV and fdh schemes, see Eqs. (35) and (36), have been 
ignored. According to the remark at the end of Section 2.4 this is correct if the A terms are 
also ignored in the real corrections. 
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4.1 Parton distribution functions in DRED 

The results Eqs. (37) and (57) have been given without taking into account that 
the subtracted partonic cross sections have to be multiphed by PDF. In dred 
it might seem natural to distinguish partonic cross sections with ^ or ^ in the 
initial state and convolute them with different PDF. We will show that this is not 
required and that we can use Eq. (57) for initial state full gluons g, convoluted 
with just one PDF even in DRED. In particular, there is no need to introduce 
unphysical PDF for finding a ^ in a hadron. 

In the strict spirit of DRED it is correct to consider independent PDF for g 
and g and write a hadronic cross section as a sum of the form 

fg/H ® c^'3"ps(^i • • •) + fg/H ® doTpsigi ■■■) + fq/H ® dcTpsiqi • • •) + fq/H O d&psiqi ■■■)■ 

(60) 
The partonic cross sections da{gi) can be constructed in the same way as Eq. (57). 
They are individually finite and satisfy 

Y^ Y^ ^DRED(ai) CJDRED(a2) j^ ,^ ^ \ J- / \ ^ r t \ 

} } 7 — T 7 — T rfa-ps(ai,a2;. ..a„) = da-ps(ai,a2;...a„). (61) 

"^^ ^-^ t^DRED(ai) C^DRED(a2) 



a\ a2 



All the PDF in Eq. (60) would be obtained by performing a fit at one particular 
factorization scale /xq and then using Altarelli-Parisi equations to evolve them to 
any other scale /i. 

The central point is that the unphysical PDF fg/H is of the order e, and hence 
its contributions to both the hadronic cross section, Eq. (60), and to the evolution 
of the other PDF are of the order e and thus negligible. 

In order to prove this we start by noting that since in other regularization 
schemes one gluon PDF fg/H is sufficient it is possible to arrange the fit in dred 
such that fg/nif^o) = at the starting scale /xq- The evolution is given by the 
Altarelli-Parisi equations, generalized to include g: 

f) f fq/Hiz) \ ^ rl^, f Pq^q Pg^q Pg^q \ f fq/uH) \ 

"V frn H ^i T ''r ^^ 'r 'rV^ • "" 

^ V fg/niz) J -Jz ^ \ p^^- p.^- p-^- y y /-/^(^) y 

where we have suppressed the fi dependence of as and fui/H and have used the 
short-hand notation Pi^j = P^J^JI^Az/^). The evolution of fg/niz) gets contribu- 
tions from Pq^g X fg/H, Pg^g X f^/H and Pg^g x fg/n- They are 0(e) x 1, C(e) x 1 
and 1 X 0{e), respectively, confirming that fg/H^jj) = 0{e) for all scales /j,. This 
in turn implies that the contribution to the evolution of fq/n and fg/n due to 
fg/H is also 0{e). The situation is in fact very similar to the contributions due 
to quarks if there were Np = e fiavours. Finally, da^^sidi, 0,2] ■■ ■ an) as appearing 
in Eqs. (60) and (61) is finite for all initial states ai,a2 separately. Hence the 
contribution of the unphysical PDF fg/n in Eq. (60) is 0{e). 
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This confirms that DRED can be used throughout all parts of calculations for 
hadronic cross sections at one loop, without the need for unphysical PDF. It is 
correct to consider only one gluon PDF fg/H and the combined cross sections 
da{ai,a2] ■ ■ .an) from Eq. (57) also in DRED. In other words, Eq. (1) is correct 
in all schemes, including dred, if the sums over all parton types include only the 
full gluon g and quarks g, q (and possibly further, massive partons). 

4.2 Summary of practical computations in DRED 

Let us finally summarize how to do a next-to-leading order calculation in practice 
in DRED or any other scheme. The main point is that only for the ultraviolet 
renormalization and the collinear counterterm a split of g into ^ -|- ^ is required. 

Virtual Corrections: To obtain the virtual corrections we start by computing 
A^Rs (fli • • • fln) with ai G {g,q,q}- For the actual calculation of the one-loop 
diagrams we do not need to split the process into many different parts as in 
Eq. (46), but can compute directly with g. The structure of the ultraviolet 
counterterms depends on the RS and the symmetries of the underlying theory. If 
DRED is used in a supersymmetric context, even the counterterms can typically be 
computed without the split, and usual multiplicative renormalization is sufficient 
to generate the counterterms. Using HV or CDR in supersymmetric theories leads 
to the complication of non-multiplicative, supersymmetry-restoring counterterms. 
In non-supersymmetric theories, determining counterterms in DRED requires the 
split g = g + g. 

Once we have the renormalized one-loop matrix element, A^rs(oi • • • On) we 
obtain rfcr^g* by integration over the phase space. Again, a split of da^g^ as in 
Eq. (61) is not required. The split g = g + g for external gluons becomes useful 
if we want to express the singularity structure or RS dependence in a simple way 
as done in the next subsection. 

Real Corrections: The real corrections in DRED can be obtained in a straight- 
forward way by directly integrating the 4-dimensional tree-level matrix elements 
containing only 4-dimensional gluons g and (anti) quarks. Likewise, in CDR we 
have to integrate the D-dimensional tree-level matrix elements. Regarding the 
real corrections in HV and FDH we remind the reader of a subtlety related to 
unitarity (see also Ref. [15]). At first sight it might appear that there is no differ- 
ence in da^g^ for schemes where the tree-level matrix elements are evaluated in 
four dimensions, e.g. between da^y^ and dal'j^}^. However, this is not correct. In 
order to maintain unitarity, in the singular regions initial and final state partons 
have to be treated in the same way as partons in a closed loop, i.e. as "inter- 
nal" . Thus in HV and FDH it is not correct to simply integrate the corresponding 
four-dimensional real tree-level matrix elements over the phase space. In partic- 
ular, Eqs. (30) and (32) contain 0{e) terms in FDH and HV and result in finite 
differences between da^^^^^ and (io'pp^j'j, even though the tree-level matrix elements 
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agree. 

In principle this procedure leads to the A terms in Eq. (35) because in the 
HV and FDH schemes incoming and outgoing splittings diflFer. As discussed in 
Section 2.4 it is possible to redefine the results by ignoring the A terms if the 
same is done in the collinear counterterms. 

Collinear Counterterm: The collinear counterterm is given by Eq. (58). It 
realizes the MS factorization scheme independent of the RS used for the compu- 
tation. For DRED, we stress that for the term given explicitly on the r.h.s. of 
Eq. (58) the partons ai . . . a„ as well as a^ are never g or g separately, but can 
always be combined to g. For the virtual parton Oj, however, it is important to 
treat g and g separately [7] . This requires the use of the splitting functions given 
in Eqs. (98) and (99). For the HV and FDH schemes, the A terms can be ignored 
in accordance with the computation of the real corrections. 

4.3 Translation rules between different schemes 

In the following we describe how the results in the various RS are related, making 
use of the split g = g + g as appropriate. We will focus on the virtual corrections; 
similar results for the real corrections and the collinear counterterms can be 
trivially obtained from Eqns. (52), (54), (58). Starting from our renormalized 
result in CDR, A1cDR(ai . . . a„), written in the form of Eq. (23), we can obtain 
■Miiv{o,i ■ ■ ■ an) simply by replacing the D dimensional (colour-linked) Born terms 
by the corresponding 4-dimensional expressions. No further change is required, 
since 7cDR(ai) = 7Hv(ai). 

The only difference between the HV and FDH scheme on the other hand does 
come from the differences of the 7^3 terms, which have been explained in Sec- 
tion 1.3 and quantitatively given in Eq. (26). For a process with ^g gluons, #g 
massless (anti) quarks and ^Q massive (anti) quarks the difference is 

M^pL{ai ■■■an)- A^Hv(ai . . . a„) (63) 

_ "^ A^(o) f^ ^ \ ^^ 41^ fljivjax) - 7FDH(a. 
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M^^^niai ■■■an) ^ #a; ( 



•^s . .(0) / 



^e{g,ij,(5} 



In the second line the influence of the different 7rs for all parton types is made 
explicit, in the third line the result in brought into a compact form. 

The difference between FDH and DRED is obtained by taking the difference 
between Eqs. (48) and (23). We can bring it into a simple form by using that we 
can write the tree-level quantities in FDH in a DRED-like form, 

Mf^n{ai ...ai...an) = ^ A^DRED(ai . . . a^ . . . a„). (64) 



24 



The difference is then governed by the factors 7dred(^) — iFonig) and 7dred(^) — 
7fdh(5')- As explained in Section 1.3 the first of these vanishes. The second is 
non-zero, 

7dred(^) - 7fdh(^) = ^ {N, - 2 TfNf) , (65) 

o 

and is present for every g in the initial or final state, see Eq. (48). Exploiting 
also that at leading order a^ = ctg, we obtain 

M'oLuiai ■■■an) - M'piijiiai . . . a„) (66) 

as iFunjg) -7dr.ed(^) -^:^ ^~/r-iA A^(o) f- -\ 
In e ^ — ' 

{a} 

as 12TfNf- N^^ .JO) , m 

= ^- ^ }_^A^UB(ai...a„)|^^_.^^. 

i=i 

Again, in the second line the infiuence of the different 7rs relevant for g is made 
explicit. In the third line we used that for processes with at least one ^, i.e. with 
#^({a}) > 1 we have A^dred ~ e; therefore we neglected the 0(e) terms from 
Eq. (65). The notation A^DHED(ai . . . a„)| _ implies that all gluons except 

gluon ij are 4-dimensional gluons. Thus in the final expression on the r.h.s. of 
Eq. (66) we sum over all processes where one 4-dimensional gluon (^ at a time is 
replaced by a ^. 

The transition rules Eqs. (63) and (66) ignore 0[e) terms and are given for 
pure QCD processes. However, they can easily be generalized to other processes, 
involving e.g. photons or massive partons, simply by using the corresponding 
explicit expressions for 7rs. 

5 Conclusions 

The main result presented in this paper is that dred can be used for the cal- 
culation of cross sections at NLO, even for processes with hadrons in the initial 
state. Problems related to factorization, as reported in the literature [8-10], can 
be avoided by taking into account the appropriate, generalized factorization in 
DRED [7]. We have shown explicitly how to use DRED together with an arbitrary 
factorization scheme. In particular, the conventional PDF [17, 18] in the MS- 
factorization scheme can be used. Also we have given explicit rules on how to 
transform separately the various parts of the hard partonic cross section, Eq. (2), 
from DRED to other RS. This completes the previously known set of transition 
rules between CDR, HV and fdh [13-15]. It is thus possible to use different RS 
for different parts of the calculation which might help simplifying the explicit 
computations. In this context we also reiterate the distinction between fdh and 
DRED. According to the definitions of the RS given in Section 1.2, at one loop 
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FDH is equivalent to the scheme DR used e.g. in Refs. [13-15] but differs from 
DRED used e.g. in Refs. [3,11,12]. 

The sahent feature of a consistent use of DRED is the spht g = g + g- In 
practice, this spht does not significantly complicate calculations. It is needed 
mainly for the correct treatment of the collinear limit of squared matrix ele- 
ments. Thus it affects the collinear counterterm and the phase-space integration 
over the singular, collinear region. The modifications regarding the former are 
shown in Eqs. (35) and (54). For the phase-space integration, the usual proce- 
dures have to be slightly modified. For the method presented in Ref. [30] for 
example, the collinear singularities in the real corrections are singled out using 
distributions, enforcing the collinear limit of the real matrix element squared. If 
this method is to be used together with DRED this simply means that the proper 
collinear limit, Eqs. (51) and (53), has to be taken. For the dipole subtraction 
method [27] additional dipoles with g are required. These can be obtained making 
minor modifications of existing dipoles, similar to the corresponding adaptation 

to FDH [32]. 

In some cases, the split g = g + g is also required for the ultraviolet counter- 
terms, since e.g. the couplings gqq and gqq renormalize differently. This seems to 
be a disadvantage of DRED. On the other hand, one of the advantages of DRED 
is that in supersymmetric theories no supersymmetry-restoring counterterms are 
required (in many practical cases; for a recent discussion see Ref. [3]). In this 
case, also couplings with g and g renormalize identically and renormalization 
is actually simpler in DRED than e.g. in CDR. This facilitates the use of the 
DR-scheme for supersymmetric parameters which is used in a wide variety of 
calculations [33]. 

Thus, DRED is a RS which is well compatible with supersymmetry and which 
can be realized with minimal modifications compared to CDR and used for an arbi- 
trary cross section at NLO. In the past, following the examples of e.g. Refs. [9,34], 
many predictions for supersymmetric processes at hadron colliders were calcu- 
lated using CDR in spite of the required supersymmetry-restoring counterterms. 
In the future, similar calculations can alternatively be carried out using DRED, 
which can lead to simplifications with the present, better understanding of DRED. 
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A Collinear limits and sum rules 

In this appendix we study the colhnear hmit of squared matrix elements and 
derive the associated sphtting functions and 7rs terms. Even though most of the 
results presented here are well known, we repeat them for the reader's convenience 
and to fix our notation and conventions. 

Following Ref. [15], we consider a slightly off-shell massless outgoing parton 
ai{pi) that splits into massless on-shell partons ak{pk) and ai{pi). The momenta 
are parametrized as 

Pi = .p.. + .l4^. (67) 

with p^ = 71^ = (k± -p) = {k± ■ n) = 0. The invariant mass of the incoming parton 
is p? = 2[pk ■ pi) = —k^/lzll — z)) and vanishes in the collinear limit A;^ -^ 0. 

To start with we consider the particularly interesting case where the parent 
parton is a gluon. Denoting by A^ the amplitude of the parent process, stripped 
of its polarization vector e'^{pi), we can write the collinear limit of the full process 
as 

M^hl ■ ■ ■ 9iPk),g{pi) . . . n + 1) ^^i^' i^^ 5_ A^ V<^X,, A- , (69) 

Sr fx Sri 

where m and n are colour labels. After averaging over /cf, the operator 'p<"^^^ 

° *= -L' ^ gi^gkgi 

is proportional to the metric, which in accordance with Eq. (10) corresponds to 
the polarization sum ^e^e*'^ . Explicitly we find 

^~' + z{l-z)] , (70) 



' g*^gg y^l " 


= -g^''{2N, 


n-r(^) = 


-- -r'(2iV, 


' g*^gg \^) 


= -r'(2iV, 


' g'^gg v^/ ~ 


= -r^'c^N, 




~r (2iv, 



1 — z z 

z 1 — z 



1 — z z 

z l-z 2 



1-z z D-2 

z l-z 2 



l-z z D 

z l-z 



z{l-z)j , (71) 

zil-z)^ , (72) 
zil - z)] (73) 



l-z z 

The interesting point is that in dred we get a combination of g'^'^ and g'^'^. Thus 
the collinear limit has to be written as a sum over two terms as in Eq. (30). In 
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the spirit of DRED, we can further disentangle the sphtting operator and write 

^< DRED flU q-y< DRED flU ^. ^< DRED flU ^. ^< DRED ^U ^. ^< DRED ^U (7A\ 

9*^99 ~ 9*^99 ^9*^99 ~^ 9*^99 ~^ 9*^99 ^ ^ 

-9'-' (2iVe) i=-^ zil -z)- r (2iV,) 1^ - r (2iVe) ^ . 

The sphtting functions P^^gg can be read off of Eqs. (70) - (74) simply by drop- 
ping the polarization sum. Performing similar calculations for all other possible 
splits we find the following results: 



n^^T, = P^^ = ^Nj-- + ^ + zil-z)) , (75) 



z \ — z 

g'^gg ^ g'-^gg ^ \ \ _ 2, ' Z 



4 
D~2 



n^'".^ = P^^ = P^^ = Tf{1- 7^^^(1 -z)) , (77) 



P^^Z = P^^ = ^^ ( T^ + ^(1 - ^) ) ' (7^ 



2z 
1~ 



n^™ = n^°r = Cf[-^-, + {1-z)] . (79) 



The results for HV are always identical to the ones for CDR, 



p< CDR 


P<HV 


LED are 


given 


by 




)< DRED 
9*^99 


= 2N^ 


4- 
D- 


-2 ^ 


)<DRED 


= 2Nc 


/l 

v 


_ ^X 


9*^99 




-)< DRED 
9''~^<19 


= Tf, 






J< DRED 

q*^99 


= Cf- 


4- 


^(1- 



^0) 



-z), (81) 

(82) 

(83) 

z). (84) 

The remaining splitting functions can be obtained by -Pjf_^|.;(-2) = Pii^iki^ ~ z). 
The splitting functions appropriate for the split of an incoming parton P^^*i{z) 
can be obtained through the crossing relation Eq. (33) with 

Cc;Hv(fi') = ^FDh(5') = ^DRED(fi') = 2 (iV^^ - 1), cjRs(g) = 2 A^c, 

uJcDR{g) = cuuREDig) = {D- 2){Nl - 1), 

Cc'DRED(^) = (4-D)(iV2-l). (85) 
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We note that in CDR and in dred there is no difference between the sphtting 
functions for incoming and outgoing partons, i.e. Eq. (34) holds with all A terms 
equal to zero. Thus, Eqs. (81) and (82) and Eqs. (83) and (84) are not inde- 
pendent. For the HV and fdh scheme the A terms do not vanish. In the fdh 
scheme, for example, the different coefficient of the three terms on the r.h.s. of 
Eq. (76) lead to A^^"^ 7^ 0. In the HV scheme, the different ujiiy{g) 7^ ^cDR(fi'), 
together with Eq. (80), are the origin of the non- vanishing A"^ terms. 

The splitting functions -Pjf_!^|; are defined only for z < 1. We define the full 
splitting functions through the relation 

Pt-^kii^) = J^—^P^^^Mi^) + S.klns{a^) 6{1 - z) , (86) 

where we made use of the standard -|-prescription. The factors 7Rs(aj) and thus 
Pi^ki ^^6 determined by requiring that the momentum sum rules 

/' dz z [P-;,(z) + 2Np P,%(^)] = 0, (87) 

/' ^^ ^ [n'"..(^) + n";.(^)] = 0' (88) 

are satisfied in all the schemes CDR, HV, fdh, i.e. also taking into account terms 
of higher-order in e if appropriate. Eq. (87) determines 7rs*(5') and Eq. (88) 
determines 7rs*(q')- For dred, the sum rules given in Eqs. (87) and (88) have to 
be generalized in an obvious way, since we also have to take into account g: 

t d^ ' K-^I(^) + PF^Is(^) + 2^^ ^^-!^?(^)] = 0, (89) 

Jo 

fdzz[P-ZT,{z) + P----{z)] = 
Jo 

f dz z [P--(.) + P--(^) + P--(^) + P^Zf^iz)] = 0, (90) 
Jo 



fdzz[P-^ 
Jo 



^) + P^^TAz) + 2Nf P^^TM] = 0. (91) 



g^gg^ ' ' g^gg^ ' ' "^'''^g^gq 



As before, Eqs. (89) and (90) determine 7dred(5') and 7dred(q') respectively, while 
Eq. (91) determines 7dred(^)- The results are 

lliV, {3D-8)TfNf {10 -D)Cf 

Icunig) = — 3(D-2) ' 7cDR(g) = ^ , (92) 

_ (6D - 13) N, _ (SD - 8) TfNf ^ ,^Cf .„,. 

7pdhI^J- 3p_2) 3(Z^-2) ' ^™«^^^- 2 ' ^ ' 

7dREd(^) = 7FDH(fl'), 7DRED(g) = 7FDH(g), (94) 

7dred(^) = 2N, - TpNp. (95) 
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and 7Hv(aj) = TcdrIoj)- Expanding these results and taking into account all terms 
to 0{e) we obtain the results given in Eq. (26). We note that these results are also 
consistent with the quark- number conservation sum rule as well as with Eqs. (31) 
and (52). In FDH and DRED they also satisfy the supersymmetric relation 

i'^^. + SiV^P;-, = P^Zl + P^Zl, (96) 

pDRED iniyr pDRED pDRED , pDRED fQT'l 

if we set N^ = Cp = 2 TfNf- 

Finally we mention that the DRED splitting function used in Eqs. (54) and 
(55) are defined as 



pDRED 
9-^99 


. _ '^DREd(^) pDRED 1 '^DREd(5') pDRED 

~ 1,) (n) ^^^^ 1,1 (n) ^^^^ ' 


(98) 


pnnEB 
9^99 


. _ '^DREolfi'j pDRED , ^DREd(^J pDRED 

" CO^nUg) '^'' CUdred(^) '^'' ■ 


(99) 


Using the explicit results above we find 




p<DRED 

9^9*9 ~ 


o.. (D-2 z l-z D-2 \ 


(100) 


-Ac^ 2 l-z ' z ' 2 ^^^ ^>) ' 


p<DRED 

9^9*9 ~ 


o.. 4-D/ z \ 


(101) 


2 \l-z J 


p<DRED 

q^g*q ~ 


^ [2{l-z) , D-2\ 


(102) 


\ z 2 J 


p<DRED 

q^g*q ~ 


r ^-^ z 


(103) 



for the splitting functions used in Eq. (53). 

B Examples 

B.l gg -^ qq 

The process gg — > qq with massless quarks has been computed long ago at one 
loop [35] in CDR and was one of the processes used to determine the relations be- 
tween the HV and FDH scheme [13]. The one-loop matrix elements M.Y<s*{gi 9] Qi O) 
were found to be related as given in Eq. (23). This and the related process with 
massive quarks was also at the centre of claims regarding problems with factor- 
ization in DRED [8-10]. The factorization issue related to the real corrections for 
these processes has been solved in Ref. [7]. Here we focus on some aspects related 
to issues with DRED, starting with the virtual corrections. 

The calculation of the one-loop diagrams is straightforward and we stress once 
more that there is no need to disentangle g into g + g in the explicit calculation 
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-NA/VW^ 



M^u^l{g,g;q,q) 



M^o^tl{9,9;q,q) 



p — ^ 



o — ^ 




\ 



— ^ 



,0 — ^ 



-6- 



Q 




Figure 6: Ultraviolet counterterm diagrams due to coupling renormalization 
-^dreMS,^;?,^) (upper line) and M\l^il{g,g;q,q) (lower line). Dark vertices 
represent counterterms 5Z°^^° and bright vertices stand for counterterms 5Z^ 



'DRED 
'9 



of the one-loop diagrams. The only issue in the computation of A^dred(5', 9, Q, q) 
is renormalization. For massless quarks we only have to consider coupling renor- 
malization. In CDR, HV and FDH this simply amounts to adding the counterterm 



f(i) 



((1) 



-MksU^, 9; q, q) = -MksU^, 9: q,q) + 2 5Zf * M'^Ug, g; q, q) , (104) 



(0) 



where Z^^* is the RS dependent coupling renormalization factor (in the MS 
scheme). In dred we have to split the counterterm contributions as 

M'hA9,9;q,q) = M'DA9,9;q,q)+M'DA9^9;q,q) (105) 
+ A^dred(^, g; q, q) + A<dr'e*d(^, g\ q, q) 

and renormalize all four parts on the r.h.s. separately. For this we need the 
coupling renormalization factors Zg^^° for the gqq coupling and Zg^^° for the 
gqq coupling. They are well known [24] and read^ 



^^DRED 



^^DRED 



asCr {-11 + e)Nc + 4:TFNF 

47r e 

a^ cr / 1 



47r e \2N, 



6 
^ + TpA^i. + e Finite J 



(106) 
(107) 



The finite 0{as) term in Eq. (106) is required because we use the MS and not 
the DR scheme. The divergent part of Eq. (107) is determined by requiring the 
cancellation of UV singularities in the off-shell gg Green function. The finite 
0{as) terms in Eq. (107) would have to be determined by a renormalization 
scheme. However, they will not affect the final result. As is to be expected, 



^We have set as — ae in these results. This is allowed since we are working at one loop. 
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this allows us to perform the calculation without specifying a renormalization 
scheme for the unphysical gluons. The counterterm M.oREr)i9,9',(l,(l) is sim- 
ply given by 2(5Z°^^° A^dhed(^,^; g, g). As illustrated in Figure 6, the coun- 
terterm M^jIA9,9;Q,<i) is given by (5Z^°^^° + 6Z°''^'')M^^Lu{9,9;q,q), while 
■MDREoi9,9',C[,Q) is not proportional to the corresponding tree-level amplitude. 
We explicitly verified that after renormalization the one-loop matrix element 

M^uLd{9, 9; Q, q) = M^uLb{9, 9; Q, q) + M'okeI{9, 9; Q, q) (108) 

in DRED is related to the other schemes as given in Eq. (66). 

The calculation of the real matrix elements is trivial. They are simply the 
four-dimensional results, Eq. (44), and the corresponding real cross section can 
be obtained in DRED by integrating these matrix elements over the phase space. 
The only remaining and main issue is the factorization of the initial state collinear 
singularities. According to our discussion, Eq. (58), it is clear that we will have 
to add 



j_.coll 

"'''dred,ms 

as_cr 
2tt e 



{9,9;(i,(i) 

dz 



(109) 
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_ pMSe ( ^ 



dcT^Luigizpi), 9{P2);q, q) 



+ 



pDRED; 

9^m ' 



z)-PY'lL[z 



+ {p, 



*qq 



*qq\ 



pMSe ( •- 
^g^qqV'^, 



da 



(0) 



iqizpi),g{p2);9,q) 



M 



d(ToREu{q{zpi),9{P2);q,9) 



M 



+ P^^fg{z) da'^'UKzpi), 9iP2y,q, q) + {1 -- 2} . 



The conversion to the MS scheme requires the terms 

TpNp 



pMSe 
9^99 



pCDR 
9~*99 



fpCDR 1 
I 9-^991 _ 



JCDR, 
9-^99^ D~>4 



6{l-z), 



pMSe 

9^m 



-- pCDR 

~ 9^m 



;iio) 
;iii) 



3 

-K-^^Jz.^4 = -^^^2.(l-.). 

It is the last term on the r.h.s. of Eq. (109) that is non-standard and deserves 
special mention since it resolves the issue regarding the seemingly non-factorizing 
corrections in DRED. 

Let us add a comment on why the factorization problem of Refs. [8-10] was 
found in the context of the process with massive quarks, discussed below, rather 
than the one with massless quarks. The reason is the fact that in the present, 
massless case, the DRED cross sections for the g and g initial states happen to be 
equal, 

rfilED(^(pi), fl'(p2); g, q) = da^^lEu{g{pi)^9{P2); q, q)- (112) 

Hence, in Eq. (109) one can combine the terms in the first and the last line to 



' DRED V 



P!^^^.^(z)da)^Lr,{g{zp,),g{p,y,q,q), 



and the process is seen to factorize even without distinguishing between g and g. 
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B.2 gg^QQ 

Problems to reconcile factorization with dred were first mentioned in the context 
of this process with massive final state quarks [8]. As explained in Ref. [7] and 
the present paper, the factorization problem disappears if g and g are treated 
as separate partons in formulas such as Eqs. (53) and (54) or in the last line of 
Eq. (109). The reason why the apparent problem has been found only in the 
massive process gg -^ QQ is not related to quark masses but to Eq. (112), which 
happens to hold in the massless case. 

With this in mind, the massive process can be treated in the same way as 
the massless one. The only additional complication in the case of massive quarks 
is to consider the RS dependence of the mass renormalization Z^ and external 
wave-function renormalization Zq^ for massive quark lines. The RS dependence 
of these renormalization factors has been considered before (see e.g. Ref. [36]) 
and, using for example the pole scheme to define the mass of the heavy quark m, 
can be summarized as follows: 

An e \f^ J V ^ 

as cr ^ fm'^\'^ f 3 



Zr = Z^^ =1 + ^^-^0,1'-^) (---4), (114) 



l + l^-^H^l (---5) , (115) 



and Zq^ = Z^. As mentioned in the main text, 7(Q) = Cp is RS indepen- 
dent [14]. We have verified by explicit calculation that using the RS dependent 
coupling renormalization and Eqs. (114) and (115) for the mass counterterms and 
wave-function renormalization, the RS dependence of the virtual corrections take 
the form as given in Eq. (46). To use dred throughout in the calculation of this 
process we simply have to use the correct RS dependent collinear counterterm as 
given in Eq. (109) and fold the hard partonic cross sections with the standard 
PDF in the MS factorization scheme. 

B.3 gg ^ h 

In this example we consider the production of a Higgs h through gluon fusion 
in a hadronic collision. While this process is relatively simple at one-loop it is 
complicated enough to illustrate all main points discussed in the main text. The 
interaction of the Higgs with gluons is given by the Lagrangian 

£j = ^(7h/itr(F^^F^,) , (116) 

where the coupling g^ has mass dimension —1. In DRED we have to distinguish 
between the coupling for ggh, denoted by g^ and the coupling for ggh, denoted 
by cjh. At tree level the two couplings are the same, but they differ at higher 
orders. 
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The only process that contributes at tree level is g{pi) g{p2) 
elements are given by 



h. The matrix 



M'~^^{g,g]h) = gluj^sig) 



'12 



;ii7) 



with g = g in CDR,HV and FDH and g G {g,g} for DRED and -MoREol^jfi'; h) = 0. 

For the calculation of the 0{as) corrections to g g ^ h the distinction between 
gh and gh at tree level will be relevant for the renormalization. Importantly, 
however, for the non-trivial part of the explicit calculation of the virtual and real 
corrections we can set gh = gh and we do not have to distinguish between g and 
g in loop diagrams. 

First we discuss the virtual corrections and how their RS dependences arise. 
The explicit calculation of the two non- vanishing one-loop diagrams in Feynman 
gauge results in the following unrenormalized one-loop matrix elements: 



M!il{gi,g2]h) 



M'^l 
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i9i,h;h) — cr 

Ztt 



where 
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0{e' 



:il8) 



;ii9) 



^v 2 

In order to obtain the counterterms we only need to perform a renormalization 
transformation of the couplings, gh -^ Zghgh, gh -^ ZghCjh- We use the MS scheme 
to define the renormalization constants in all RS. The results read 

asCr -llN^ + ATpNp 



'7CDR 



/VDRED 

^DRED 
^gh 



^gh 



^gh 



1 + 



A-n e 

«£_££ 

47r e 



-II + e)Nc + ATFNF 



1 + 



«£_££ 

47r e 



-4:Nc + 2TFNF + eYimie) . 



(120) 

(121) 
(122) 



As expected, gh renormalizes like the square of the strong coupling, and the 
difference between Z'^f^^ and Z^^^^ is in agreement with the corresponding scheme 
difference of a, [13,37], see also Ref. [24]. Writing Z^^ = 1 



5Z^^ we have 



M^^lig, g; h) = M^il{g, g;h) + 2 6Z^t -^rsU^, g; h) 



gh 

for CDR, HV and FDH. In the case of DRED we obtain 



,W 



M\,Lu{g,g]h) = M 



(1) 



grM 



I •'*'Idred 



{g,g;h)+2 6Z^ 
{g,g;h)+2SZ^rM 



(0) 

DRED 

(0) 

DRED 



i9,9;h) 
{9,9;h) 



(123) 
(124) 



Neglecting terms of 0{e) this entails 



iW 



M\^>{g,g;h) 



a.cr (0) 



2% 



M'^i{g,g;h) 



2Nr 



£12 




(125) 
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with AJ'4 = A;;^* = 0, AZ\^ = NJ3 and A^^*^^ = 2NfTf/3. Thus, after adding 
the counterterms, these expressions are in agreement with the general formula 
Eq. (46), with the finite RS independent part given by A4iii{g,g; h) = g^ {N^ — 
1) Np 8^2/6- The scheme dependences also exemplify the formulas discussed in 
Section 4.3. 

Turning to the calculation of the real corrections, we also have to take into 
account the processes with (anti) quarks in the initial state. The corresponding 
matrix elements in CDR are given by 

M^U9,9;h,g)= (126) 



,2,1 AT f T^T'i 1\ \ ^ 



(D-2)(Si2-Sl4)^(Sl2-S24; 



2 



glAT^as N,{Nt - 1) > v- ^^v^^^ -^-^/ v^^^ ^^^^ _ 4, 



S12S24S14 



12 



cycl 

M^^Ug,q;h,q)= (127) 

gf^AnasiFiJyc - 1) 7^ , 

M'§U<l,Q;h,g)= (128) 

glAna^TpiN^-l] 



^2 ^,{D- 2){sl^ + 4) + 2{D - 4)s24Si4 



2s 



12 



where the sum in Eq. (126) is over all cyclic permutations {pi ^ P2 ^ ^^4} 
and the corresponding matrix elements in HV, FDH and DRED can be obtained 
by setting D ^ 4. Eqs. (127) and (128) are related by crossing. 

We have now all the necessary matrix elements at hand to compute cross 
sections. In this example we restrict ourselves to the total cross section a. The 
cancellation of singularities and the RS dependence is the same for any infrared- 
finite observable. 

Defining x = M^/su = p'i/su we write the leading-order cross section as 

a^h99)^^^h9,9;h) = —S{l-x){M^^^{g,g-h)) . (129) 

S12 

Note that in DRED we have <TDRED(fl'fl') = o-dred(^5') = o-DRED(^fi'), analogously to 
Eq. (112). This is another example of the special case discussed at the end of 
Ref. [7]. Obviously, the virtual corrections can be obtained by the same formula 
if we replace A^rs by A^ks . Defining the function 

-"-(-^'©'^(-^(t?)"-!-^^^!' U30) 

the explicit results read 

<s\99) = —{M^nk9, 9; h)) [S^-"' + 5{l - x)A-f ) . (131) 

S12 
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To obtain the real corrections we parametrize the phase space with the help of 
the variable y defined such that 

si4 = ^(1 -x){l-y); S24 = ^(1 - x)(l + y). (132) 

The real corrections can then be written as 

.real/ ^ X _ (ISvr)"" ^\l-2e [ j„,f, „,2\-e/ i\AO) , 



<f(«i«2) = ^T7Y3^^r2(l-^) y ^dy{l-y')-^{M^J{a,,a2;h,as)) (133) 

and evaluated using distribution identities like 

1 _ 1 ^/i ^A , 1 o /ln(l-a;) ^ ^ ^^^ ^ 



-6{l-x) + - _-2e ^ '-] +0{e 



(l-a;)i+2^ 2e ' ' (l-a:)+ \ 1-x 

= -^S{l-x) + I+{x)-2eL+{x) + 0{e^). (134) 

In CDR and DRED the integration is straightforward, while in HV and fdh the 
distinction between internal and external gluons as defined in Section 1.2 leads to 
a complication. As discussed there, in HV and fdh soft or collinear gluons have 
to be treated not as g but as g and g, respectively. In order to ensure this correct 
treatment, one can subtract the collinear limit of the integrand in Eq. (133) and 
integrate it separately, either explicitly or using e.g. the dipole formalism as in 
Ref. [15]. Alternatively, one can split the integrand using distribution identities 
for (1 + y)~^~'^ similar to Eq. (134) and treat each term as appropriate. 
Ultimately, defining the functions 

S^^-\gg) = AT, (^^^^ ' (^1 5(1 _ :,) _ 1 (/4^) _ :,(2 - x + x^)) 



+ y (x - 1)^ + 8 L4x) (x^ - a; + 1)' - y (5(1 - x) j 
- 2L+{x) (1 - x) (2x - 2 - x^) J , 



(135) 



(136) 



the real cross section can be written as 

a 



(gg) = ^{M'ii{g,g;h)){S-^\gg) + A^J!,) , (137) 

Sl2 



<s\9Q) = —{M^nk9,9;h)){S-%q) + AI^,), (138) 

With A'aU = AZ = AZ.^n = 0, Af„H = -4Ar,x2(l - x), AgL, = A^ = and 
ApDH = Ad'^ed = —Cpx'^. The matrix element Mrg^s{q,q]h,g), Eq. (128), will 
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not produce any singularities upon integration over phase space. Thus cr'^s^lqq) 
is finite. 

The remaining ingredients needed for the hard partonic cross section are the 
colhnear counterterms suitable for the MS factorization schemes. They are con- 
structed according to Eq. (58) and read 



da^"^^ 



Ms,iis\99 ) 
as 
2tt e 



(139) 
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pMS € I 






d(T^:Hs{9q) 
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^'^ Idz 



pRS / 
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z)jda^^^{g{pi),g{zp2);h) 
z)da^^^{g{pi),g{zp2);h) 



where 



pMSe 
<1-'9Q 



pCDR 



fpCDR 1 



—eCp z , 



(141) 



in addition to Eq. (110) is taking into account the conversion to the MS factoriza- 
tion scheme. In order to present the explicit results for the colhnear counterterms 
we introduce the functions 



S'^'^'igg) 

S^°\gq) 



-4N, ixH+{x) + (1 - x){l + x^)) + -/5o5(l 



X] 



-Cf[2-2x + x^) . 



(142) 
(143) 



The colhnear counterterms can then be written as 



,coll 



0"r 



{99) 



a 



^{M'^i{g, g; h)) {S^^'\gg) - A|| - A™* 5{l - x)) 
S12 



<^{9q) 



«, 



S12 



f{0) 



{M'::i{g,g;h)){S^^\gq)-^^^,). 



It is now easy to see that the subtracted, partonic cross sections 

a{gg) ^ [<y^"\gg) + <yZ\99) + <ylf{99) + <{99)\n. 
H9q) = [<s\9q) + ^^\9q)]n^, , 

Hqq) = [<s\qq)]n-., , 



(144) 
(145) 



(146) 
(147) 
(148) 



are finite and RS independent. In the sum all the RS dependent Apg terms drop 
out, and in all RS we obtain the well-known result that can be found e.g. in 
Ref. [38]. 



37 



To summarize the calculation in dred: Both the one-loop diagrams and 
the real corrections and phase space integrals can be computed in a straight- 
forward way, using only full gluons g. The split g = g + g has to be used in 
the evaluation of the collinear counterterms rfcr^°" and in the computation of 
the UV counterterms. In the end, the cross sections d'{gg),a{gq),a{qq), as well 
as o-{qg),d-{gq),a-{qg),cr{qq) are RS independent. The hadronic cross section is 
obtained by convoluting them with the standard parton distribution functions 
obtained in the MS factorization scheme. 
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